Fast adiabatic-like spin manipulation in a two-electron double quantum dot 
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We apply the transitionless quantum driving method to control the electron spin of a two-electron double 
quantum dot with spin-orbit coupling by time-dependent electric fields. The x and y components of applied 
electric fields in each dot are designed to achieve fast adiabatic-like passage in the nanosecond timescale. To 
simplify the setup, we can further transform Hamiltonian by z axis to design an alternative speed-up adiabatic 
passage, without using the applied electric field in y direction. 
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Introduction — Spintronics aimes at fast and robust 
spin control in nanostructurei~— . There are several meth- 
ods to manipulate spin accurately, such as electron spin 
resonance induced by magnetic field oscillating at the 
Zeeman transition frequency^ and electric control with 
spin-orbit (SO) coupling^. Recently proposed techniques 
of "shortcut to adiabaticity"— ~— motivate us to achieve a 
high-fidelity control in a single quantum dot 9 (QD) and 
a two-electron double QD^S. In a single QD, we applied 
inverse engineering method to design a fast and robust 
protocol of spin flip in the nanosecond timescale, based 
on the Lewis-Riesenfeld theory^. Furthermore, in a two- 
electron QD, more freedoms of the applied electric fields 
provide the flexibility to control the spin from one ar- 
bitrary state to the target state, by controllable Lewis- 
Riesenfeld phases. A different shortcut is provided by 
tansitionless quantum driving^, reformulated by Berry 
and equivalent to counter-diabatic control proposed by 
Dcmirplak and Rico^. This technique was originally uti- 
lized to control the spin in the fast adiabatic-like way 
by the applied magnetic fields^^. Shortly afterwards, 
it has been extensively applied to various quantum sys- 
tems like two-level or three-level atomsi, Bose-Einstein 
condensates in accelerated optical lattices^ and electron 
spin of a single nitrogen- vacancy center in diamond^. 

In this Letter, we use transitionless quantum driving 
to design the external electric fields for rapid spin ma- 
nipulation in a two-electron double QD in the presence 
of a static magnetic field and spin orbit coupling. Differ- 
ent from Berry's calculation^, we will apply the electric 
fields and take advantage of spin-orbit coupling, since the 
time-dependent electric fields are easy to be generated on 
the nanoscale by adding local electrodes^. In addition, in 
a single QD, there are only two controllable parameters, 
that is, x and y components of the electric fields, so that 
it is difficult or even impossible to produce the desired all- 
electrical interaction by transitionless quantum driving^. 
However, transitionless quantum driving is applicable in 
a two-electron double QD, as there exist more freedoms 
with four controllable parameters, x and y components 
of the external electric fields for each electron in double 
QD. To simplify the experimental setup and reduce the 
device-dependent noise, we can further apply the concept 
of multiple Schrodinger pictures^ to find an alternative 
shortcut with only x component of the applied electric 
fields in our system, by z rotation of Hamiltonian. 
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FIG. 1. (color online.) Schematic diagram of a two-electron 
double quantum dot in the presence of the external electric 
fields. 



Hamiltonian — Two electrons are confined in a dou- 
ble QD, described as a quartic potential, where they are 
isolated by Coulomb blockade, illustrated in Fig. [TJ The 
spin-dependent Hamiltonian is fltotai = H s + H- m t, where 



H s — JS\ ■ S2 
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(i) 



Here j = 1,2 represent the electrons 1 and the electron 
2. Sj — (Tj/2, are the spin operators of two electrons. 
The Zeeman term Aj = g/isBj, where /ib = h\e\/(2moc) 
is the Bohr magneton, g is the Lande factor with nega- 
tive value (g < 0), like in GaAs and InAs, and mo is 
the mass of the free electron. B\ and B2 are the static 
magnetic fields applied to the electron 1 and 2 in z direc- 
tion, respectively. For simplicity, we set B\ = B 2 = B, 
so that A = Ai = A2 = gjisB. In the presence of the 
applied magnetic field, the lowest four eigenstates of the 
system can be expressed by singlet and triplet for S = 
and S 1 = 1 in the basis of \S,S Z ). If the energy difference 
between the singlet |0, 0) and the lowest one of the triplet 
|1, 1) is much less than the gap between the singlet and 
the triplet J, which means | J + A| <§; J, we focus on the 
state transition between these lowest two states |0, 0) and 
|1,1), as shown in Fig. [1] By choosing 1 , 0) = (1,0) T 
and |1, 1) = (0, 1) T , we can write the reduced Hamilto- 
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nian H s in the form of 2 x 2 matrix 
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-J 







II- (2) 

-J + A 

The interactions between the electric field and the elec- 
tron H m t are expressed as, 



(3) 



where the vector potential Aj (t) are related to the exter- 
nal electric fields, £\{t) = —(l/c)dA\/dt and £ 2 {t) = 
— (l/c)dA2/dt, and vj are the spin-dependent veloc- 
ity operators. We consider the SO coupling including 
structure-related Rashba (a) term and bulk-originated 
Dresselhaus (/3) term for [110] growth axis, 



(4) 



so that the spin-dependent velocity operators become 

i ]=/3<T*-ao>j, (5) 



[H soc ,y 3 } = aa] 



(6) 



Therefore, the total spin-dependent Hamiltonian -fftotai 



tl total — ^ 
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To rewrite the symmetric Hamiltonian, we shift one 
quantity Z Q = -J/4 + A/2 - (e/3/2c)(Af + A%), and 
finally obtain Hamiltonian -fftotai = ZqI + H, where 



H 



Z X + iY 
X-iY -Z 



(12) 



with Z = (1/h) [- J - A + (e0/c)(Af + Afj\. The solu- 
tion to the Schrodinger equation of H differentiates from 
that of -f/totai by the factor exp[— i J Zo(t')dt'], while the 
energy level |0, 0) and |1, 1) are shifted by Z$. The states 
after the shifting arc denoted by |1) and | — 1), respec- 
tively, and their populations remain unchanged as the 
ones of the previous states, |0, 0) and |1, 1). 

Transitionless fast spin tranfer — Our aim is to trans- 
fer the spin from | — 1) to |1) totally during a reasonably 



short time duration. The form of Hamiltonian H in Eq. 
[T2l tells us that Y and Z are the functions of A\ and 
A% and X is the function of A\ and A\. Different from 
the Hamiltonian of one electron confined in a single dofeS,, 
the transitionless quantum driving can be applicable to 
the spin control in a two-electron double QD, as we can 
figure out how the Hamiltonian H (including reference 
Hamiltonian Hq and counter-diabatic term H 1 ) is imple- 
mented by corresponding electric fields. We may take the 
reference Hamiltonian H as 



Hn = - 



Z iY 
-iY -Z 



(13) 



driven by Af and A% . The example of a double QD 
of GaAs-based structure is considered below, where g = 
—0.44 and the static magnetic fields are B\ = B 2 = 3.43 
T. The energy gap between the singlet and the triplet is 
J = 0.1 mcV, so that | J + A\/J = 0.12 < 1 with the 
above parameters. 

With the help of reference Hamiltonian ((T3)) . we can 
write down the instantaneous eigenstates, |x±), satisfy- 
ing Hq\x±) = E±\x±), where the instantaneous eigen- 
values are E± = ±h\f Z 2 + Y 2 /2, and the instantaneous 
eigenstates are 
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, (14) 



with the mixing angle 6 — arccos[Z/ (Y 2 + Z 2 )] and ip 
7r/2. Once the adiabaticity condition^ 



ZY - YZ 



{Y 2 



< 1 



(15) 



is fulfilled, the state 1^°}, the solution to the Schrodinger 
equation of Hq, evolves from |\E ,0 (0)) = |x±(0)) and fol- 
lows the adiabatic approximation 



= exp 



dt'E±(t') 



\x±(t)). (16) 



Otherwise, transitions between |x±(0)) will occur. To 
implement population inversion, from | — 1) to |1), 
along one of instantaneous eigenstate, \\+(t)), we set 
the ansatz of the vector potential A* = Aotanh[(t — 
cijtf)/(wjtf)], where etj, Wj describe the change rate of 
A j and j = 1,2. To fulfill the initial and final states, 
Y(0) = Y(tf) = should be fixed, which means at the 
initial and final times Af and A^ are equal to each other. 
Meanwhile, the mixing angle 9 goes from ir to 0, crossing 
the point tt/2 during the interval (0, t/). With this strat- 
egy, we produce the reference electric fields £ J , displayed 
in Fig. [21(c). For the following comparison, we first show 
the dynamics of populations for the instantaneous eigen- 
states, Pi n (t) = |(l|x+(t))| 2 andi^t) = |(-l|x+W>| 2 
(seen in Fig. [Ha)). In practice, this process is not adi- 
abatic, so the populations of exact solution, x E ,0 (t), of 
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FIG. 2. (color online.) (a) Time evolution of the populations 
P{ n (t) (solid blue line) and P_\(t) (dashed red line) as the in- 
stantaneous eigenstates of Ho, which coincide with the pop- 
ulations Pi(t) and P-i(t) as the solution to the Schrodinger 
equation of H. (b) Time evolution of the populations P\{t) 
(solid blue line) and Pl 1 (t) (dashed red line) as the solution 
to the Schrodinger equation of Ho, showing that this is not an 
adiabatic process, (c) The applied electric fields in the x di- 
rection £f (solid blue line) and £% (dashed red line), and the 
additional two electric fields in the y direction with the differ- 
ence £ = £% — £% (dot-dashed orange line) drive the popula- 
tion inversion. Other Parameters: tr = 2 ns, ha = 1.2 x 10 -6 



meV 



h/3 = 0.3 x 1(T 6 meV 



Hamiltonian H are obtained as P°{tf) = |(l|*°(t/))| 2 = 
0.76 and P°i(t f ) = \{- l\^°(t f )}\ 2 = 0.24, (seen in Fig. 
[2](b) ) , which are not consistent with that of the instanta- 
neous eigenstates. Of course, the adiabatic passage can 
be realized by extending tf and increasing the electric 
fields, respectively. For example, if we prolong tf = 14 
ns and keep the previous A*, the process will become 
adiabatic, and P\{tf) = 0.9999 is finally achieved. On 
the other hand, Pi(tf) — 0.9999 can be also achieved, 
when the magnitude of £J are increased by 11.5 V/ m 
and keep tf = 2 ns. 

Next, transitionlcss quantum driving will provide sup- 



plementary time-dependent interactions Hi that can- 
cel the diabatic couplings of a reference process H , 
and make the reference process fast and adiabatic- 
like. The supplementary counter-diabatic term Hi = 
E±\dtX±)(x±\ 



Hi 



X 
X 



(17) 



driven by A\ and A\, where X = 6 = (YZ - YZ)/(Z 2 + 
Y 2 ). As a result, the solution <I>(t) to the Schrodinger 
equation of H = Ho + Hi becomes exactly the adia- 
batic approximation of Hq. The corresponding dynam- 
ics of the populations, Pi(t) = |(l|*(t))| 2 and P-i(t) = 
|( - l|*(i))| 2 . The populations Pi(t) and P 2 (t) coincide 
with P{ n (t) and -P™ W respectively, as shown in Fig. HJa). 
The difference between two additional y components is 
A V D = A\ — A\ = HO/ {s/2ea) and the corresponding time- 
dependent function of £ V D = £\ — £f is plotted in Fig. 
[2jc). The maximal magnitudes of £ V D is 0.94 V/m. Ob- 
viously, they are much less than the increasing values 
in magnitude of to achieve the adiabatic process, as 
mentioned above. This implies that the transitionless 
quantum driving can really speed up the adiabatic pro- 
cess. As a matter of fact, £ V D , as the function of 8, is 
related to Y and Z. The shorter time is, the larger value 
of £ V D is required. To implement £ V D easily in the exper- 
iment, we need the smooth function of £ V D , therefore in 
general £f and £f should not vary very dramatically. 

z-axis rotation — In reality, the electron spin is sub- 
ject to the device-dependent noise, which could be the 
amplitude noise of the electric fields^. It can be quite im- 
portant, especially when the electric fields are relatively 
weak. From the above analysis, we find that four con- 
trollable parameters, £J and £j, x and y components of 
the electric fields for each electron in double QD should 
be applied. If y component of the electric fields can 
be reduced, we can decrease decoherent effects resulting 
from the device-dependent noise. To this end, we can 
apply the concept of multiple Schrodinger pictures, and 
make unitary transformation of Hamiltonian H by z-axis 
rotation^. We write down the dynamical Hamiltonian 
as follows 



H= H - 
2 



Z it 

-iQe-^-^/V 



J e »(0-7T/2) 

-Z 



(18) 



where tan<^> = Y/X and Q = V X 2 + Y 2 . By applying 
the unitary transformation^ 

h ( e *(?W2) n \ 
U =2{ e-^)J' 

which amounts to a rotation around z axis by the angle 
7r/2 — <f>, we calculate the new Hamiltonian H' = U^(H — 
K)U with K = ihl/W , and finally obtain 



H' 



h ( z 



iQ 



-iQ 



(20) 
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FIG. 3. (color online.) (a) The populations P™(i) and 
P-i(t) as the solution to the Schrodinger equation H' . (b) 
Comparisons between £f n (solid blue line) and Sf (dashed 
red line), (c) Comparisons between £f" (solid blue line) and 
£f (dashed red line). Other parameters are the same as in 
Fig.H 



without a x term. We should notice that the dynamics 
of Hamiltonian H = Hq + Hi and H' is not the same 
(the populations are the same because of z rotation). 
However, the Hamiltonian H' is equal to the original one 
H = Hq + H i at t = and t / , which guarantees that the 
initial (final) states of H and H' coincide. So the Hamil- 
tonian H' can provide an alternative way to implement 
the shortcuts to adiabaticity. According to the Hamilto- 
nian H' Eq. (|20[) . we may acquire two new controllable 
parameters, £f n and x component of the electric 

fields, since Z + (£> and Q are the functions of the sum and 
the difference of Af ™ and " , respectively. The solution, 
^ n (t), of the Schrodinger equation of H' can be solved 
numerically, and the populations Pi(t) = |(l|\I/ n (i))| 2 
and P'Mt) = l( - il*"^))! 2 : ar e shown in Fig. E| (a). 



At the final time, -P™(£/) = 1 and the population is com- 
pletely inverted. The new electric fields only in x di- 
rection are shown in Fig. |3] (b-c) with some corrections 
compared with the previous ones £f and £f • 

Conclusion — We propose the shortcuts to manipulate 
the spin states formed in a two-election double QD by us- 
ing transitionless quantum driving. The Hamiltonian H 
is divided into two parts, the reference process Hq, driven 
by Af and A\ , and the supplementary time-dependent 
interaction H\ , driven by A\ and A\ . By applying x and 
y components of electric fields for each electron, the spin 
system follows exactly the adiabatic approximation of the 
reference Hamiltonian Hq, in the time scale of nanosec- 
ond. In order to simplify the setup, and decrease the 
device-dependent noise effect, we further transform the 
Hamiltonian by z axis and obtain the new Hamiltonian 
implemented only by x component of electric fields. This 
provides an alternative shortcut to realize the fast and 
adiabatic-like spin control. We hope these results may 
lead to the applications in spintronics and quantum in- 
formation processing with the state-of-the-art technique. 
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